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1. INTRODUCTION 

Research by Thukaram et al. [1] forward and backward sweep power flow is used. According to 
Lohia et al. [2] genetic algorithm is intended for optimum setting and sizing of the capacitors for distributed 
system. The objective function used is cost minimization due to real power losses as well as reactive power 
compensation in system. Supply systems is normally self-possessed of feeders, different transmission 
erections, are circled. The best significant moments of expending radiated feeders are that many patrons are 
affected by the disaster of a particular element. Reliability valuation and optimization are two important 
modules to performance overall lessons for supply networks [3], [4]. In the past distribution systems (DS) 
were given lower importance in reliability studies compared to other networks. Currently, huge quantities of 
developments regarding dependable and reasonable operation of supply systems are being settled in [5]-[7]. 
Numerous practices constructed on systematic and typical practices have been advanced for the consistency 
evaluation in supply systems [8], [9]. 

Nowadays penetration of the distributed generation (DG) with photovoltaic (PV) on DS was 
maximized drastically [9]-[15]. When DG is installed it can solve numerous commercial and ecological 
problems, the reliability of the system also improves [13]. If the DG is installed at unsuitable location may 
prime to surge in power losses, and even the voltage magnitude will be affected [12]. The suitability of 
supply networks by reliability indices average failure rate A, average outage period r and yearly outage period 
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U. The reliability are capacity indices of essential prominence in dependability study. Power system 
dependability acts a significant role in the DS that is more than eighty percent of customer disturbances 
occurs due to disasters in the DS [16]. Research by Gupta et al. [17] presents a performance on optimal 
reliability strategy of electrical DSs. DG connection in DS, the constancy and safeguard the distribution 
network performed a detailed exposition in [18]. 

Junior et al. [19] explained that a new method was proposed to regulate the interval outcomes of a 
power system load flow for three-phase unbalanced distribution networks depends on three-phase current 
injection method (TPCIM) on state variables were deemed at rectangular coordinates. At presented method, 
the active and reactive powers at every busbar were modeled as interval values to account for its inherent 
uncertainties and Krawczyk operator was used on the power flow equations to deliver reliable interval three- 
phase outcomes. In addition, the usage of interval extensions and angular rotation system to overcome the 
overestimation issues related with interval solutions was presented. Chen et al. [20] elucidated a three-stage 
relaxation-weightsum-correction (TSRWC) based probabilistic reactive power optimization system to 
contract with the uncertainty and correlation of doubly fed induction generator (DFIG) fed wind generators 
on distribution network. Primarily, the inverse Nataf transformation, an enhanced point estimate system was 
presented to transform the uncertain wind speeds. Ali et al. [21] elucidated a novel system for optimal DS 
planning with DG systems. The penetration of DG in power systems is a method that has numerous benefits, 
like peak savings, reduction of losses and improvement of voltage profile. It is also intended to maximize 
system reliability, stability, and security. The main objective of optimal distributed generation (ODG) was a 
assurance to attain the aids mentioned above to maximize the general efficiency of the system. In the study, 
an improved wild horse optimization (IWHO) algorithm was presented as a novel metaheuristic system to 
solve optimization problems on electric power systems. 

Numerous optimization systems have been proposed in the optimum location of the reactive power 
compensation devices. The proposed method is executed in the IEEE-33 radial distribution feeder. The 
efficiency of proposed method is established in MATLAB environment. 


2. RADIAL STABILITY INDEX 

Voltage stability index was proposed [22] depending on the existence of solutions for real as well as 
reactive power at the getting end. Relative strength index (RSI) in [23] is given in (1). Figure 1 describes the 
electrical equivalent. 


RSI(r) = 2V2V2 — Vj — 2V,2(PR + QX) — |ZI2(P? +Q?) í) 


RSI(r). voltage constancy index of node r (r=2,3,4 ...... ) 

For steady process of the radial distribution networks, RSI(r) > 0 for r =2,3,4....... 
Where, 

Vs: implies sending end voltage 

Vr: implies the getting end voltage 

Z: implies the impedance of the branch 

R: implies the resistance 

X: implies the reactance 

P: implies the real power load 

Q: implies the reactive power load 


P+jQ 


Figure 1. Electrical equivalent 
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2.1. Reliability indices 


The system average interruptions frequency index (SAIFI): 


SAIFI = ZEN (2) 
N; 


l 


Where A; is implies failure rate of it” load and Ni represents the count of customers supplied at load. 
The system average interruption duration of index (SAIDD, 


SAIDI = SLN (3) 


i 


Where U; indicates unavailability of load I that is product of failure rate and repair rate of i” load. Customer 
average interruption duration index (CAIDI) is, 


SAIDI 


CAIDI = (4) 
SAIFI 
Average energy not supplied index (AENS) is expressed by, 
—_ Total Energy not Supplied by the system 

AENS = Total number of customers served (5) 
2.2. Derivation of proposed delta analysis system 

A linearized power flow is expressed by, 

AP] _/H_ Ny [A@ 

k = (u L) lav (6) 

AP = J,,A0 + Ji2AV (7) 

AQ = J240 + J22AV (8) 
Rewriting (6) we get, 

A0 = JītAP — Jita AV (9) 
Substituting 40 from (9) in (8) we get, 

AQ = Joy AP = Jair Ji2AV + Jo2AV Uiii — J22)AV = hjit A- AQ 

AV = Voit 12 = J22 Uzt AP — AQ] (10) 
Then, ARSI (r) is determined by, 

ARSI(r) = [4V, — 4V? — 4V, (PR + QX) [Davie — Jez) Vali AP — AQ]] (11) 


In (11) gives the relationship between the alteration in RSI index with variation in P and Q. In (11) is utilized 
to detect the sensitivity of voltage stability indices using the reactive power injection at a specific bus. 


2.3. Steps for optimal location of reactive power compensation devices 


In order to maintain the voltage profile at each bus within an acceptable range, improve voltage 


stability, reduce power losses in lines, increase system reliability and security, network optimal usage is 
crucial with the modernization of power grids. These can be accomplished by integrating reactive power 
compensation devices into distribution or transmission networks. The steps for optimal location of reactive 
power compensation devices is described as follows, 


Run the power flow to obtain sending end and receiving end voltages. 

Calculate the radial stability index given in (1). 

From (11) obtain the matrix which gives the relationship between ARSI(r) and AQ (i.e.AP = 0). 
Consider a limit to the radial stability index. The buses which have lower index value then the limit is 
considered as most prone to instability. 

Obtain the reduced sensitivity matrix by keeping only buses, which is more prone to the voltage 
instability. 
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— Route LP optimization technique using MATLAB to get location and size of reactive power 
compensation device. 


2.4. Algorithm for optimal location of distributed energy resources 
Small generating plants known as DG are connected to consumers through DSs to enhance voltage 
proficiency, voltage profile, stability, decrease of power losses, and economic benefits. The foregoing 
advantages are attained by optimally placing DGs. The steps for optimal location of distributed energy 
resources (DERs) is described as follows, 
— Run the power flow to obtain sending end and receiving end voltages. 
— Calculate the radial stability index given in (1). 
— From (11) obtain the matrix which gives the relationship between ARSI(r) and AP (i.e. AQ = 0). 
— Consider a limit to the radial stability index. The buses which have lower index value then the limit is 
considered as most prone to instability. 
— Obtain the reduced sensitivity matrix by keeping only buses, which is more prone to voltage instability. 
— Route LP optimization technique using MATLAB to get the DERs location and size. 


3. RESULTS AND DISCUSSION 
3.1. Test case 1: optimal location of reactive power compensation devices 

The method derived in section 1.1 is executed on IEEE-33 bus system. Table | displays voltages 
profile of IEEE 33 bus. The bus voltages before reactive power compensation are ranked and according to the 
ranking, the bus which is ranked -1 is said to be under voltage collapse and the other buses are ranked 
accordingly. Reactive power compensation devices are injected for voltage stability and we can see the 
voltage is improved. 


Table 1. Comparison of voltage profiles before and after placement of reactive power compensator (RPC) 
Bus No. Voltages before RPC Ranking Voltages after RPC 


1 1 33 1 
2 0.995972408 32 0.997263756 
3 0.976792583 27 0.985102731 
4 0.966543652 25 0.980166008 
5 0.956398113 22 0.975670554 
6 0.931134752 21 0.97275448 
7 0.926300556 19 0.980862179 
8 0.907576634 16 0.986661427 
9 0.898875669 15 0.993279639 
10 0.890800666 12 1.001054439 
11 0.889606272 11 1.001262548 
12 0.887522748 10 1.001622252 
13 0.879019895, 6 1.013062311 
14 0.875850947 5 1.020361776 
15 0.873884726 4 1.024511889 
16 0.871979941 3 1.027425034 
17 0.869155953 2 1.034629845 
18 0.8683 10362 1 1.035521801 
19 0.995270827 31 0.996563093 
20 0.9905 18897 30 0.991817401 
21 0.989583078 29 0.990882811 
22 0.988736298 28 0.990037145 
23 0.971993533 26 0.980345454 
24 0.963062486 24 0.971492315 
25 0.958608723 23 0.967077563 
26 0.928496683 20 0.971488216 
27 0.924989344 18 0.969895293 
28 0.909338255, 17 0.966288 128 
29 0.898092062 14 0.964126355 
30 0.893222104 13 0.962666983 
31 0.8875 19848 9 0.96600412 
32 0.886265327 8 0.96701465 
33 0.885876585 7 0.968242216 


The system comprises of 33 nodules, 32 sections, and 5 tie switches. The immoral arrangement 
takes exposed lines 33 to 37. Minimal active and reactive load of system is 3.715 MW, 2.300 MVAr [24]. 
Figure 2 shows the bus voltages before and after RPC. The optimal location and size of reactive power 
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compensation devices is provided to the Table 2. Power loss before and after placement of RPC is given on 
Table 3. 


E Voltages before RPC m Voltages after RPC 


1,1 


1 
9 
8 
7 
6 
5 
4 
3 
2 
1 
0 


123 45 6 7 8 9 101112131415 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 
BUS NUMBERS 


0 
0 
0 
0 
0 
0 
0 
0 
0 


BUS VOLTAGES 


Figure 2. Bus voltages before and after RPC 


Table 2. Location and size of the reactive SVC 
Size of the reactive 


Sl.No. Bus No. 


power device in KVAr 

1. 9 0 
2, 10 0 
3. 11 243.8 
4. 12 350 
5. 13 350 
6. 14 350 
Ti 15 350 
8. 16 350 
9. 17 350 
10 18 350 
11 29 0 
12 30 350 
13 31 350 
14 32 350 
15 33 350 


Table 3. Comparison of losses before and after SVC 
Active power loss __ Reactive power loss 
Before RPC 295.92 231.68 
After RPC 200.84 163.36 


From Table 3, it can be observed that out of 15 no buses which are below the index limit assumed 
only the placement of RPC at 10 buses only. From Table 4, it is perceived that the real power loss decreases 
by 69% and reactive power loss by 70.5%. According to Lohia et al. [2], after placement of capacitors active 
power loss decreases by 66.51%, reactive power loss decreases by 60.81%. Even though the losses due to 
delta analysis are little towards higher side, proposed number of reactive power devices are less than the 
method proposed in reference [2]. So, the expenditure incurred on installing the devices will come down. 
Delta analysis method is a very easy to implement as it uses conventional radial power flow and linear 
programming technique only. It may be viewed as Figure 3 that there is a substantial improvement in the bus 
voltages after reactive power compensation. 
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Table 4. Comparison of voltages profile before and after placement of DERs 


Bus No Voltages before Voltages after Bus No Voltages before Voltages after 
` placement of DERs placement of DERs ` placement of DERs placement of DERs 

1 1 1 18 0.86831036 0.9591487 

2 0.995972408 0.997293509 19 0.99527083 0.9965929 

3 0.976792583 0.985186788 20 0.9905189 0.9918473 
4 0.966543652 0.980167951 21 0.98958308 0.9909128 

5 0.956398113 0.975474003 22 0.9887363 0.9900671 

6 0.931134752 0.962207633 23 0.97199353 0.9804299 

7 0.926300556 0.958882046 24 0.96306249 0.9715776 

8 0.907576634 0.952291977 25 0.95860872 0.9671632 

9 0.898875669 0.950780069 26 0.92849668 0.9610455 
10 0.890800666 0.950960136 27 0.92498934 0.9596009 
11 0.889606272 0.950112144 28 0.90933826 0.951816 
12 0.887522748 0.950581342 29 0.89809206 0.95065576 
13 0.879019895 0.952128798 30 0.8932221 0.95053653 
14 0.875850947 0.952675298 31 0.88751985 0.95051215 
15 0.873884726 0.954653347 32 0.88626533 0.95050306 
16 0.871979941 0.956777269 33 0.88587658 0.95052623 
17 0.869155953 0.958653112 


3.2. Test case 2: optimal location of distributed energy resources 

The method derived in section 1.1 is tested on IEEE-33 bus system. Table 4 indicates the bus 
voltages before and after placement of DERs. Optimal position and size of the DERs is given in Table 5. 
Power loss previous and next placement of reactive power compensation devices is given in Table 6. It can 
be observed from Table 6 that even though 15 no of buses exceed the limit, optimal location of DERs is only 
8 numbers. Active power loss and reactive power losses decreases by 43.50%, 43.2% which are lower 
compared to losses mentioned in reference [2]. It concludes that the optimal location of DERs would reduce 
the system losses and improve the bus voltages. 


Comparison of voltages before and after placement of DERs Table 5 Location and size of the DERs 


j j | | | | | | | | | Size of the reactive power 
m 


© 
© 


Sl.No. Bus No. device in KW 


| | i 9 0 
0.7 Before DERs 2. 10 0 
| E After DERs | 3. 11 0 
se! | 4. 12 0 
= 0.5 | | 5. 13 0 
@ 0.4 6. 14 0 
ae | | J: 15 141 

| | 8. 16 200 

0.2 | | 9. 17 200 

0.1 10 18 200 

| | 11 29 0 

K 5 10 15 20 25 30 12. 30 167 

Bus No. 13. 31 200 

14 32 200 

Figure 3. Bus voltages before and after placement DERs 15. 33 200 


Table 6. Comparison of losses before and after DER 
Active power loss Reactive power loss 
Before DER 295.92 200.84 
After DER 128.73 86.843 


4. RELIABILITY ASSESSMENT 

Once after finding the optimum position and dimension of DG utilizing delta analysis method, the 
reliability system with DG is simulated by using MATLAB site. To study efficiency of proposed system, a line 
contingency is created by removing lines between buses 8 to 9 and between buses 28 to 29. Following the above 
said line contingencies the buses from 9 to 18 and buses 29 to 33 will be fed by the DG sets already installed. 
Reliability indices given in the section 1 are calculated before contingency and after contingency. 

From the Table 7 it can be observed that the indices values have decreased which indicates that the 
reliability of system is improved. Average root mean square (ARMS) and relative errors of output calculation 
shown on Table 8. Efficiency comparison is provided on Figure 4. 
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Table 7. Comparison of voltages profile before and after placement of RPC 
Index Before DG After DG 
SAIFI 0.384242424 0.281818182 
SAIDI 0.29530303 0.212727273 
CAIDI 0.768533123 0.75483871 
AENS 31065.15152 _23793.93939 


Table 8. ARMS and relative errors of output calculation 


i . Average Maximum 
Error terms (%) Parameters [25] Proposed [25] Proposed 
+ Mean calculation 0.08 0.071 0.013 0.07 
SD calculation 1.06 1.01 1.39 1.111 


ARMS error Cumulative distribution function (CDF) estimation 0.085 0.075 0.124 0.108 


Efficiency (%) 


S Base 

GA 
® Proposed 
1 PSO 


Figure 4. Efficiency comparison 


5. CONCLUSION 

In this paper, the presentation of delta analysis method is utilized for enhancing voltage stability and 
moderate power losses in system by using optimization technique, the location and size of reactive power 
compensation devices and DG units are used for improving voltage stability. Also, the study is mainly 
concerned with assessment of reliability under contingency condition by placement of DGs. Thus, the 
objectives of this paper are realized. It is detected from the delta analysis method is prevailing and easy tool 
to optimally locate both reactive power compensation devices and DGs. 
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